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Abstract 

We combine the NRQCD formalism for the inclusive color singlet and octet pro¬ 
duction of charmonium states with the parton and the ACCMM model, respec¬ 
tively, and calculate the momentum distribution in the decay B —> J/ip + X. Ne¬ 
glecting the kinematics of soft gluon radiation, we find that the motion of the h 
quark in the bound state can account, to a large extent, for the observed spectrum. 
The parton model gives a satisfactory presentation of the data, provided that the 
heavy quark momentum distribution is taken to be soft. To be explicit, we obtain 
Ep = 0(0.008 — 0.012) for the parameter of the Peterson et al. distribution function. 
The ACCMM model can account for the data more accurately. The preferred Fermi 
momentum pp = 0(0.57 GeV) is in good agreement with recent studies of the heavy 
quark’s kinetic energy. 


*This work was supported in part by the Bundesministerium fiir Bildung, Wissenschaft, Forschung und 
Technologie (BMBF), 057DO93P(7), Bonn, FRG, and by the US Department of Energy under contract 
DOE/ER/01545. 


1 






1 Introduction 


The ARGUS and CLEO collaborations reported resnlts on inclnsive B meson decays to 
J/iIj, where they identified a sizable component of decays with three or more particles 
in the final state |l|. In a recent pnblication p[ CLEO presented an analysis based on 
their sample of data which is an order of magnitnde larger than those of previons stndies, 
corresponding to a redaction of errors by a factor of 2.4. As a result the group finds the 
direct branching ratio B{B J/ijj + X) = (0.80 ± 0.08)%. 

In the context of the color singlet (wave function) model this process is related to the 
decay of a 6 quark in the B meson, at short distances, into a color singlet cc pair plus 
other quarks and gluons. The c and c quarks generated by this method have almost equal 
momenta and reside in the ^Si angular momentum state. Several authors performed the¬ 
oretical studies of the branching ratio B{b —>■ J/iIj + X), to leading order in using the 
color singlet approach Q. The measurements indicate that this approach underestimates 
the data by roughly a factor of three. The inconsistency remains when the next-to-leading 
order perturbative corrections are included ii- Thus it is interesting to consider gener¬ 
alizations of the color singlet model. 

The nonrelativistic QCD (NRQCD) factorization formalism developed by BBL P] 
makes possible a systematic treatment of soft gluon effects in the inclusive heavy quarko- 
nium production. It allows for the creation, at short distances, of a heavy quark and 
anti-quark pair in a color octet configuration which subsequently evolves into a physical 
bound state through the emission and absorption of soft gluons. 

Color octet contributions to inclusive charmonium production have been recently in¬ 
vestigated in hadronic reactions at collider and fixed target experiments 0, I, I, 0, in 
e’''e“ annihilations |^, 7 -nucleon reactions at fixed target and HERA energies [^, |1^, p!4[ , 
in decays at LEP |0, and in lepton-nucleon reactions [0. A comprehensive review 


of recent developments in this field can be found in Ref. [0. In our analysis, we shall 
refer to previous studies of inclusive B meson decays to J/V’ [0, 0, 0]- These articles 
considered the decay of free h quarks, including also the leading color octet contributions 
to the branching ratio B{b J/ip + X). However, the branching ratio is very sensitive to 
the numerical values chosen for the Wilson coefficients and the NRQCD matrix elements. 
In this situation it is desirable to study additional predictions of the theory. A second com¬ 
parison deals with the J/ip momentum distribution. Adopting the NRQCD formalism to 
leading order in the nonrelativistic expansion, we shall argue that the observed momentum 
distribution can be attributed, to a large extent, to the B meson bound state corrections. 

So far there is no detailed theoretical fit of the momentum spectrum. Palmer and Stech 
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have made a first attempt using a simple wave function formalism [^|. In this paper we 
investigate two different approaches and compare them with the experimental results. As 
stated above, a sizable component in the decay consists of nonresonant multi-particle hnal 
states. Consequently, for a wide range of the phase space an inclusive description based 
on quark-hadron duality is appropriate. This approach was applied extensively to the 
inclusive semileptonic decays of B mesons. 

Over the last few years, considerable progress has been made in the calculation of 
nonperturbative corrections to the lepton energy spectrum of the decay B —> ev + X 
using the operator product expansion (OPE) pT[. The latter approach involves a series 
in powers of 1/(1 — y)mi, with y being the normalized lepton energy and incorporates the 
formalism of the Heavy Quark Effective Theory (HQET) [^. However, an analysis of the 
decay B J/'ip + X using the HQET is of limited validity, since due to the smaller energy 
release the convergence of the OPE is slower than for the semileptonic decays. Therefore 
we shall resort to a formalism in which the momentum distribution of the heavy quark in 
the B meson, the latter representing the dominant bound state effect, is modeled in terms 
of parameters obtained from experiment. The motion of the h quark will be referred to as 
the smearing of its momentum or the Fermi motion within the B meson. 


Several effects may contribute to the momentum spectrum of the 

• the motion of the b quark in the bound state, 

• the shifting of the cc momentum due to perturbative QCD corrections to the b —>■ ccs 
matrix element, carried out beyond the leading logarithm approximation (LLA), and 

• the emission or absorption of soft gluons at the hadronization stage. 

In this work we focus on the dependence of the spectrum on the smearing of the b quark 
momentum in the B meson using the LLA for the short-distance matrix element. As 
far as soft gluon radiation is concerned, we account for the color flow using the NRQCD 
formalism, but neglect the momentum flow carried by the soft gluons. This approximation 
corresponds to keeping the leading nonvanishing orders in the n^-expansion for the color 
singlet and octet states, respectively (n being the relative velocity of the c and c). One 
might note, however, that the n^-expansion fails to converge near the boundaries of the 
phase space. In this region, which is sensitive to the mass difference between the partonic 
cc and the hadronic J/ip hnal state, the momenta of soft gluons become important and 
inhuence the spectrum 
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By analyzing the dependence of the J/V’ momentum distribution on the initial bound 
state corrections, we explore one of the two relevant nonperturbative effects numerically, 
keeping in mind that the soft gluons, which are radiated within the fragmentation process, 
must be included in the future. To be explicit, we investigate to what extent the b quark 
motion in the B meson can account for the observed spectrum in B ^ J/'ip + X decays. 

The paper is organized as follows. In Section ^ we outline the basic ideas of the NRQCD 
approach as they apply to the inclusive charmonium production in b decays. We adopt the 
NRQCD to the decay of a 6 —> J/'ip + X, where we emphasize that the b quark is considered 
to be free but the cc pairs from both color singlet and octet states convert into the hadronic 
J/'ip hnal state. The latter step is treated in the leading nonvanishing orders of the v^- 
expansion for the singlet and octet states, respectively. In order to account for the measured 
J/'ip momentum spectrum, we shall combine this formalism with the B meson bound state 
corrections, i.e., we consider the smearing of the b quark momentum in the meson. This 
we do in two models. In Section we present the parton model which gives a satisfactory 
presentation of the data. Then we calculate in Section ^ the J/ip momentum spectrum 
in the framework of the ACCMM model [^, which is in better agreement with the data. 
In this context, we take the Fermi motion parameter pp within the range obtained from 
theoretical studies of the b quark’s kinetic energy. Our analysis conhrms that, by including 
the leading color octet contributions to the charmonium production, one can account for 
the measured branching ratio. In addition, we hnd that the b quark motion gives a good 
representation of the spectrum. The summary and the discussion of our results can be 
found in Section |^. 


2 b —> J/'i/j X in the NRQCD Formalism 

The nonrelativistic-QCD (NRQCD) formalism, developed by BBL [Q, represents a compre¬ 
hensive theoretical framework for the study of inclusive charmonium production. Denoting 
by V the relative velocity of the two heavy constituents inside the (QQ) bound state, the 
NRQCD approach introduces the complete structure of the quarkonium Fock space to a 
given order in v'^ and thus allows for a consistent factorization of long- and short-distance 
effects. Consequently, the charmonium production rate is represented by a sum of prod¬ 
ucts, each of which consists of a short-distance coefficient, associated with the creation 
of a heavy quark and anti-quark pair in a specihc angular and color conhguration, and 
a nonperturbative NRQCD matrix element which parameterizes the 

subsequent evolution of the intermediate cc(^^~^^Lj)a state into a physical charmonium 
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bound state H (plus light hadrons). The latter step is generated through the emission and 
absorption of soft gluons. 

The separation of the distance scales (see Fig. 1), and its association to the small value 
of the constituents’ relative velocity (u^ ~ 0.23 — 0.30 for charmonium), makes possible 
the calculation of the heavy quarkonium production rate in terms of a double series in 
powers of and (numerically as(m^) ~ u^). This way the production through an 
intermediate color octet state, although down by powers of in the Fock state hierarchy, 
can be numerically relevant in view of possible short-distance enhancements. The authors 
of Refs. |T^, [^, pointed out that, for inclusive B decays to S-wave charmonia, the 
color octet production mechanism is indeed important, even though it is of order with 
respect to the basic color singlet one, because the Wilson coefficient of the color octet term 
is strongly enhanced in comparison to that of the color singlet transition. Since we shall 
use formulas for the decay rate F(6 —>■ J/'ip + X) in our analysis of the ACCMM model, 
we briefly outline the basic ideas of the NRQCD approach. 


At the b mass scale, the Hamiltonian which induces the effective b —>■ ccqf (/ 
transitions reads 







C.)c-f^Lcqf^^Lb + {C+ + C.)c-f^LT^cqn^^LT% 


s, d) 


( 1 ) 


L = 1 — 75 , where operators arising from penguin and box diagrams have been neglected. 
The renormalization (Wilson) coefficients C±{^) have been computed up to the next-to- 
leading order corrections in Ref. [^. 

In the leading logarithm approximation (LLA), the decay amplitude refers to the tree 
level matrix element of the effective Hamiltonian. Consequently, the renormalization scale 
dependence contained in the Wilson coefficient functions, C±{^)^ cannot be cancelled by 
this matrix element. While this work was being completed, a study of higher order per¬ 
turbative corrections to the matrix element oib ^ J/'ip + X decays appeared P|, in which 
a double series expansion in and the small ratio of the (LL) color singlet to octet 
Wilson coefficients was performed. Using the color singlet approximation for charmonium 
production, the authors obtained a result which carries only weak dependence on the renor¬ 
malization scale. The predicted branching ratio, B{B —>■ J/ip + X) = 0.9)*)^'^ x 10“^ lies 
well below the experimental value of (0.80 ± 0.08)% reported by the CLEO collaboration. 
This conhrms the supposition that a nonperturbative effect in the B J/ip + X transi¬ 
tion is required, which goes beyond the color singlet approximation, namely the color octet 
(NRQCD) mechanism for charmonium production. In future studies the leading order re¬ 
sult for the color singlet contribution can be replaced by the expression given in Ref. [^. In 
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this article we keep the leading logarithm approximation, assuming that the modification 
of the J/'ip momentum distribution will be small. 


In order to obtain the nonrelativistic interaction from the effective Hamiltonian of Eq. (|l|), 
one expresses the Dirac bilinears of the c and c in terms of (two-component) Pauli spinors 
^ and T], the constituents’ relative three-momentum q, and the Pauli matrices cTj. This 
allows one to consider separately states with specific angular momentum and color (for 
details see Ref. 0). To linear order in n = |q|/rUc, the short-distance amplitude A (using 
a notation in which the spinors do not carry color indices) can be rewritten as 


G 

A{a,T;c,d,eJ;q) = 


^mik, 


xsYLb [2mcA^A - r]r , 


( 2 ) 


c, d, e, / being color indices, a and r the individual spins of the charm quarks. P is the 
total four-momentum of the c and c in the laboratory frame, is the Lorentz boost matrix 
that takes a three-vector from the cc restframe to the laboratory frame. 

The decay rate to a cc pair, calculated in full perturbative QCD, is obtained from 


P(6 


cc+qj) 


1 


d^q d^Pf 

(27r)3 (27r)32Po 


{27ry6^^\pb-pf-P)\A\‘^m^^ + ... , ( 3 ) 


where the dots denote higher order terms in the v^-expansion. The integrations represent 
the three-body phase space of the final state. Using perturbative NRQCD, the corre¬ 
sponding decay rate to a cc pair of a specific angular and color quantum number state 
[n]a = (^®+^Tj)a reads 


P(6 ^ cc[n]a + Qf) = 


d^q 

(^ 


Sa[n 


(0|Of[n]|0) 


+ ... 




(4) 


Replacing the heavy quark field operators, which appear in the cc matrix elements, in 
terms of the Pauli spinors ^ and p permits a matching of the two perturbative results in 
Eqs. (^ and (|[). This determines the short-distance coefficients ^^[n], which equally apply 
to the inclusive decays to charmonium bound states (for more details see Refs. P, 3)- 
Consequently, the rate of the transition b —*• J/tp + X in the NRQCD factorization 
formalism reads 


T{b^J/'iP + X) 


ln]a 


( 0 |Oi/^[n]| 0 ) 


m 


2 

c 


1 


d^p d^Pf 
(27r)32Po (2vr)32E/ 


{2'kY 6^‘^'^ {pi 


Pf - P)^c ^ 


(5) 
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( 6 ) 


^ (0|0“[n]|0) ' ' “ ^ ' 


Note that the integration over the relative momentum in Eq. (Q) properly relates the 
short-distance 1 —> 3 particle transition to the long-distance effective 1 —2 particle one. 

Evaluating the color octet intermediate states in Eq. (H) up to order relative to the 
color singlet ‘baseline’, one obtains in the restframe of the b meson (with ~ 0 and 
\Vcs\^+\Vcd\^^l) 11111 , 11 , 

. 2 \ 2 r 


r(6 ^ J/i) + X) = ^^\Va,\^mcml 


mi 


, . I 

a ( 1 H- 7^ ) -|- b 

mi 


where we defined the NRQCD coefficients 

(0|0f/^(35i)|0) 


a = {2C+-C.] 


+{C+ + C.f 


3ml 

(0|0s'-'’'ys,)|0) (0|O,"^*(=‘Pi)|0) 

~r A 


2mi 


mz 


b = {C+ + C.) 


Mo\oZ'*{^s„)\o} 


2ml 


(7) 


( 8 ) 


(9) 


We shall apply Eqs. 0-0 in Section 4 in the context of the ACCMM model in order to 
determine the J/'ijj momentum distribution from B meson decays. 


In the case of the parton model, which we shall introduce in Section 3, we have to relate the 
short- and long-distance matrix elements squared rather than the decay rates. To this end, 
let us define the matrix element M. of the effective interaction b ^ J/'ip + X/m which the 
b ^ Qf transition is assumed to be independent of the inclusive charmonium production 
process, through 

r(. . ( 10 ) 

If we identify the total four-momentum P of the c and c with the momentum of the J/'ip 
in the final state, = P, then comparing Eqs. (|[) and (|T0|) we obtain 

\M\Hb ^ J/i. + X) = m-' 5^ ld!!h-^A!T^(0|O„"«[n||0). (11) 

[il\a 

An explicit calculation yields 

|AfP = Gl\V,,\%A{^f+mf)Ym + m,)YL ), (12) 
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where the J/V’ tensor structure reads 


= \Vc,\V,- 


“1“ (fl “1“ j^2 i 


(13) 


with a and b dehned in Eqs. (H) and (P). The momentum identihcation, described above, is 
in accordance with the BBL formalism for quarkonium production. To leading order in the 
n^-expansion, the soft gluons, radiated in the nonperturbative transition cc —>■ J/ip + X 
(and materializing into light hadrons), are assigned no energy or momentum; i.e., the 
energy dependence of the long-distance interaction is neglected, and thus the NRQCD 
matrix elements appear as universal constants. 

In Section 3 we shall apply Eq. ([l^) in the context of the parton model. For comparison 
with data we identify 2mc with the mass Mp of the J/'0, because to leading order in the 
n^-expansion = 4m^ (however, one might keep in mind that the expansion fails to 
converge near the endpoint of the spectrum). 


3 B —> J/X m. the Parton Model (PM) 


3.1 Calculation of the Differential Branching Ratio 


As we mentioned above, in the context of the NRQCD approach the nonperturbative 
charmonium production process cc —>• J/'^-1-X is assumed to be independent of the 
transition. Thus, within the leading logarithm approximation for the short-distance matrix 
element, the B J/'ip + X amplitude can be evaluated using a (generalized) factorization 
approach; i.e., the amplitude is related to a product of matrix elements of current operators 


M{B^J/i; + Xf) = 


-(2C+ - C.){X}\qf^,Lb\B){ij + X^lc^LclO) 


+ iC+ + C.){Xj\qf^^LT%\B){^lJ + X«|c7^LT“c|0) 


(14) 


with Xj + X^ = Xj + X^ = X'j. Note that, contrary to the color singlet wave function 
model, Eq. (pN]) includes the color octet transitions. 

The standard NRQCD approach outlined in Section 2 considers charmonium production 
from a free b quark decay. As it becomes obvious from Eq. (^, the corresponding transition 
b J /if + X refers to a 1 —> 2 particle phase space, because we neglect the soft gluon 
momenta. In this approximation, the J/if momentum distribution from B meson decays 








arises from the initial bound state corrections (and the momentum smearing from the 
Lorentz boost to the laboratory frame). We incorporate the corrections, attributed to the 
interaction between heavy and light quark in the initial meson, in the parton model along 
the lines developed for semileptonic B decays ^7_ . 


If we use Eq. ( 0 . which is valid at the level of tree calculations, in the restframe of 
the B we can write the decay rate as 


dT{B J/'iP + Xf) = 


Gl\V, 


cb\ 




d^k. 


■y 


' 2K, 


(15) 


{2tiYMb 

with dehned in Eq. (13). Introducing the light-cone dominance, as in Ref. [^, allows 
us to relate the transition B ^ Xf to the distribution function f{x) of the heavy quark 
momentum, 

= 4:{S^pu\-ief,p„x) f dxf{x)PB{xPB-k^ye[{xPB-k^)o]S[{xPB-k^f-mj], (16) 
Jo 

with 


SppuX 9fj,p9uX 9iiy9pX “1“ 9px9vp and ^( 2 ^) 


(17) 


- 1 - 1 , x > 0 
— 1, a; < 0 . 

The dependence of the distribution function on the single scaling variable x is a consequence 
of the light-cone dominance, since in this framework the structure function /(x) is obtained 
as the Fourier transform of the reduced bilocal matrix element between the hadronic states 

0 . 




•■B 


The Lorentz invariant width of the decay reads 

Eb- dV = ^ / dxf{x)e[{xPB - [{xPb - k^f - m]] 

X _ ^ ^ •X 


( 18 ) 


( 19 ) 


f=s,d 


X 

where 

Evaluating Eq 
momentum spectrum, 

1 dr 


a - b)PB{xPB - ky + -j^{PBky{xPBk^ - M^){a + h) 


dJk^ 

2E^i, 


^ Gl 


7271 


VMVcfl^Ml 


( 20 ) 


in the restframe of the B meson, we arrive at the formula for the J/y 

kyp 


(R ^ J/y + X) = tb 


Ts d|k^| 

X I [f{x+) + f{x_)] 


Gl 

727tMb 




E, 




( 21 ) 


{a - h) + {a + &)ypjY 


M, 




+ [fix+) - fix-)] (a + b) 


2E^ I 'ktp I 
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within the kinematical range 


M 


^/j — ^'ip ^ 


Ml , n/r2 _2 




+ 


2M, 


B 


Here we defined 


x± = 


± 


+ m 


( 22 ) 


(23) 


and adopted mj ~ 0 and | 


Mb 

P + iKdP — 1- The kinematical range for a;_ belongs to a final 
state qnark with negative energy. Therefore the corresponding terms can be associated, 
formally, with quark pair-creation in the B meson whereas the dominant terms proportional 
to /(x+) reflect the direct decay. 

To compare our result with data from CLEO a Lorentz boost has in addition to be 
performed from the restframe of the B meson to a moving B produced at the T(4S') 
resonance (|pn| = 0.34 GeV). The explicit form of the boost integral is presented in 
Section ^T| in the context of the ACCMM model (see Eq. (^)). 


We note that the dependence of the decay spectrum on the structure function appears in a 
factorized form. Thus the shape of the spectrum in Eq. (pT]) is governed by the functional 
form of the heavy quark momentum distribution. This is different for the semileptonic 
decay, in which an integral over the structure function is involved. In the following section 
we compare the predictions of Eq. (0) with the existing experimental data. 


3.2 Analysis and Numerical Evaluation 


In order to compute the theoretical momentum distribution we must fix the parameters 
appearing in Eq. (|2l|) . In the numerical analysis we use Gp = 1.1664 x 10“® GeV“^, 
Mb = 5.279 GeV, = 3.097 GeV and C+{mb) = 0.868, C-{mb) = 1.329. The (leading 
log) values of the Wilson coefficients correspond to the central value of the pole mass, 
rrib = 4.7 GeV Moreover we adopt tb = 1-61 ps for the B meson lifetime 

and \Vcb\ = 0.040. The latter value stems from a GLEO II analysis 


of the inclusive 

semileptonic B decays [| and is in accordance with the world average from inclusive and 


exclusive measurements, \Vcb\ = 0.0381 ± 0.0021 


As we mentioned above, the numerical values of the nonperturbative NRQGD matrix 
elements (0|Oa^^(^'^+^Lj)|0) = must be determined phenomenologically. 


^The analysis was based on a modified version of the ISGW model and on ACCMM Q, re¬ 
spectively, using Pp = 0.3 CeV for the Fermi momentum parameter. The corresponding uncertainty, 
among other things due to the ambiguity in the choice of we postpone to the discussion of the overall 
uncertainty associated with the normalization of the B —> Jji!) + X decay spectrum. 
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For the rest of this article, we use 


(Of*CSo)> 


1.1 GeV^, 
0.0066 GeV^, 
0.04 GeV^, 
-0.003 GeV^. 


(24) 


In addition we have to impose the (approximate) heavy quark spin symmetry to relate the 
matrix element (Og^^(^Po)) fo that arising from the i^Pi) intermediate state relevant for 
B decays, 

{oi'*epj)) ^ (2J+ i){op*ep„)}. (25) 

The color singlet matrix element Si)) has been determined from the leptonic width 

of the J/'ij) [Q, the corresponding color octet matrix element (Og^^(^S'i)) stems from the 
analysis of J/V’ hadroproduction at GDF |^. The remaining parameters (Og^^(^S'o)) and 
(Og^^(^Po)) have been extracted from a £t to data on J/'^ leptoproduction (the latter 
involving only a small dependence on the value chosen for m^. 

Note that the renormalized matrix element i^Po)) is not restricted to positive 

values, a characteristic due to the subtraction of large power divergences from the (positive 
dehnite) bare matrix element (for details see Ref. [16||). 

The numerical values of {OI^^Q-Sq)) and (Og^^(^Po)) quoted in Eq. (P^) are consistent 
with results from photoproduction , pion-nucleon reactions and, roughly, with cal¬ 
culations of J/V’ production in hadronic collisions [§. Yet one has to note that the various 
investigations of charmonium production processes within the NRQGD approach involve 
leading order theoretical expressions. The corresponding uncertainties in the determina¬ 
tion of the matrix elements set up a sizable range in the parameter space (for a detailed 
discussion see Ref. |]^). However, as we discuss next, an important feature of our analysis 
is that the momentum dependence is largely free of these uncertainties. 

In both models, the parton and the AGGMM model, the shapes of the spectra originat¬ 
ing from various intermediate states are essentially identical. This statement is 

obvious in case of the (^Ti) vector states, which involve identical Lorentz structures, and 
thus the relevant NRQGD matrix elements appear in a single constant (see Eq. (H)); in case 
of the (^Rq) scalar state the equivalence follows from the explicit calculation. Gonsequently, 
the shape of the entire J/ijj momentum distribution in B meson decays is practically inde¬ 
pendent of the exact values of the nonperturbative matrix elements. The same statement 
holds for the numerical values of the parameters occurring in the short-distance coefh- 
cients (in particular errors have to be attributed to the Wilson coefficients, due to scale 
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uncertainties, and to the value chosen for the CKM element Vcb)- 

The set of parameters mentioned above solely determines the normalization of the 
J/ijj momentum spectrum. Within our approach the shape is attributed to the momentum 
distribution of the b quark in the initial B meson. Further investigation of soft gluon effects 
is clearly needed. The hrst step we take here is to ask: can a successful reproduction of 
the data at hand be achieved without them. 


In the parton model, in the absence of direct measurements of the distribution function we 
use a one-parameter Ansatz and £x the distribution parameter by comparing our results 
with data. Referring to theoretical studies which pointed out that the distribution and 
fragmentation function of heavy quarks peak at large values of x [^, we assume, as a 
working hypothesis, that the functional form of both is similar. The latter is known from 


experiment and we shall use the Peterson functional form |36 


f{x) = W 


a;(l — x)" 


[(1 - xY + EpXY 


(26) 


with Ep being the free parameter and W the corresponding normalization constant. This 
function has already been applied in the semileptonic decays of the B meson . 


In Fig. 2a we show the distribution function f{x) for various values of Ep. The kine- 
matical range for the two arguments x± appearing in Eq. (pID reads 


M,. 


M. 


< < 1 , 


B 


Ml 

Ml 


< x_ < 






(27) 


B 


Note that the variable X- only occurs with values at which /(x_) is small. Therefore the 
corresponding contribution to the differential branching ratio is small. 

We use Eq. (^6[) to fit the measured momentum spectrum of the J/Y which was pre¬ 
sented by the CLEO group 0]. As pointed out above, the shape of the theoretical spectrum 
is determined by the distribution function, i.e., by the parameter Ep. 

A general feature of the analysis is the difhculty for reproducing the data over the 
entire range of phase space. Confronted with this problem, we lay greater emphasis on the 
appropriate description of the low momentum range (|k^| < 1.4 GeV). Within this region 
the J/Y spectrum obtains a sizable contribution from decay channels containing three or 
more particles in the hnal state, whereas the high momentum range is mainly determined 
by the exclusive two-body decays B —>■ Therefore incoherence, as a necessary 

ingredient of the PM, is fulhlled in the former region. In addition, the u^-expansion fails 
to converge in the endpoint domain of the spectrum. 
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As a result we find that if we apply small values for the parameter Sp < 0.004 — 0.006 
we cannot account for any part of the low momentum region which is underestimated, 
whereas the high momentum range is overestimated. The situation improves if we apply 
larger values for Sp = 0(0.008 — 0.012), corresponding to a soft b quark distribution (see 
Fig. 3). For Sp = 0.008 we obtain a hrst satisfactory £t for the present data. Yet there is 
still a moderate systematic underestimate. A further increase of the distribution parameter 
(with the remaining parameters kept unchanged) does not imply significant modifications 
within the low momentum range. For comparison, one might note that in the studies of 
semileptonic B decays £p was taken between 0.003 and 0.009 |2^ . 


Considering the uncertainties associated with the determination of the NRQCD matrix 
elements and, consequently, with the normalization of the spectrum, in Fig. 4 we show the 
J/'ip momentum distribution for Sp = 0.012 and a modihed value Po))/ml = —0.001 

GeVT The £t is somewhat better. However, the color octet intermediate states, with 
reasonable values of the matrix elements, improve the £t of the calculated curves to the 
observed spectrum. For comparison, and in order to point out the importance of the color 
octet charmonium production mechanism, we also show in Fig. 4 the spectrum arising from 
the color singlet (wave function) model = 0). The latter underestimates 

the data by roughly a factor of three. 

As a general result we conclude that a soft b quark momentum distribution can, to a 
good extent, account for the observed spectrum in i? —> J/ip + X decays. We shall see 
in Section [4.2| than an analogous statement holds for the Fermi momentum distribution in 
the ACCMM model. 

Finally, the errors in the data, although substantially improved, are still significant 
and a better test will be possible, when the error bars will be further reduced. With 
better data on should also include higher order nonperturbative corrections to the quarko- 
nium production, and study the influence of soft gluon radiation within the fragmentation 
process. 


4 B 


J/ip + X in the ACCMM Model 


4.1 Calculation of the Differential Branching Ratio 

A second approach which allows us to analyze the momentum distribution of J/ip in the 
inclusive decay of the B meson is given by the ACCMM model In this model the 

bound state corrections to the free b quark decay are incorporated by attributing to the 
spectator quark a Fermi motion within the meson. The momentum spectrum of the J/tp 
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is then obtained by folding the Fermi motion with the spectrum from the b quark decay. 
In Ref. the shape of the J/il) momentum distribution resulting from Fermi momentum 
smearing has been given in the context of the color singlet wave function model (and 
without consideration of the light spectator quark mass). We shall extend the analysis 
by including the leading color octet contribution to the charmonium production in the 
framework of the NRQCD factorization formalism and by comparing our results with 
experimental data. 

The spectator quark is handled as an on-shell particle with dehnite mass rrtsp and 
momentum |p| = p. Consequently, the b quark is considered to be off-shell with a virtual 
mass W given in the restframe of the B meson by energy-momentum conservation as 




(28) 


Altarelli et al. introduced in the model a Gaussian probability distribution (j){p) for the 
spectator (and thus for the heavy quark) momentum, 


0 (p) 



exp {-pVpI) , 


(29) 


normalized according to 

dpp‘^(j){p) = 1. (30) 

The Gaussian width pp is treated as a free parameter which has to be determined by 
experiment. 

One main difference between the parton model and AGGMM is that in the latter one 
must consider a b quark in flight. We therefore start from the momentum spectrum of 
the J/ijj resulting from the decay b —^ J/tjj + Xf (/ = s, d) of & b quark of mass W and 
momentum p which is given by 



dVb 

d\k^ 


(|kpl, V) 


% 


-1 


k^l\p) - \k^-{p)\ 


(b). 


e Ik 


ip\ 


|y‘’(p)|) - 9 (|k*| - kf\p)J 


(31) 


Here we have defined 


9{x) 


1 , X > 0 

0, a: < 0 . 


(32) 


Fq is the width of the analogous decay in the restframe of the heavy quark as obtained from 
the NRQGD factorization formalism. According to Eqs. (0)-(||) it reads (again adopting 
m/ ~ 0 and |WsP + iWdP - 1) 


ro = 


Iddvr 





a 




(33) 
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where a and b contain the nonperturbative effects in the charmonium production process. 

In Eq. (ED give the limits of the momentum range which results from the Lorentz 
boost from the restframe of the b quark to a frame where the b has a nonvanishing mo¬ 
mentum p, 

h*’(p) = -j^(SA±p£o), (34) 

with 

io = ^ - Mj) . £„ = + Ml, (35) 

and 7 ^^ being the corresponding Lorentz factor, 

7 ,“' = !^, E,= ^W^+p^. ( 36 ) 


To calculate the momentum spectrum of the J/ip from the inclusive decay of the B meson 
one has to fold the heavy quark momentum probability distribution with the spectrum of 
Eq. ([3T|) resulting from the b quark subprocess. Performing this, we finally arrive at the 
expression for the differential branching ratio for a B meson in flight. 


1 dV 

Ts d|k^ 


{B ^ J/iP + X) = 


fc+(|kv,|) 



|fc-(|kv.l)l 


dlK 


Pma 


KiKi) - \k.{K,\)\ 


dp p^cppp) 


dPft 

dm 


(Ikil, p) 


( 37 ) 


Here Pmax is the maximum kinematically allowed value of the quark momentum p, i.e., 
that which makes W in Eq. equal to W = (for comparison with data we identify 
2mc with the mass of the J/ip), 


Pmax = [{Ml + m% - Mlf - ^ . (38) 

The first integration in Eq. (^) results from the transformation from the spectrum for a 
B meson at rest to the spectrum for a B meson in flight, where 

m(mi) = (Es|k^| ± Ipsim) , k+{\k^\) = mm{k+{\k^\), kmax} , (39) 

Mb 

kmax being the maximum value of the J/ip momentum from the decay B —> J/ip + X in 
the restframe of the B, 


k — 


{Ml + Ml- - 4M|M2 


Srn.in. 


sp ' 


(40) 


2Mb l 

(Note that in the PM Smin = 0 for vanishing masses of the final state quark qj.) 

In the following Section we make use of Eq. (|37D to compare the model predictions with 
experimental data. 
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4.2 Analysis and Numerical Evaluation 


Both models, the parton model as well as ACCMM incorporate the bound state structure 
of the B meson by postulating a momentum and, consequently, a mass distribution for 
the heavy quark. In Section we pointed out that in the PM the dependence of the 
J/ijj momentum spectrum on the exact value of the Peterson parameter Sp (within the low 
momentum range |k^| < 1.4 GeV) is only moderate. We shall see that, contrary to this, in 
the ACCMM model the dependence of the spectrum on the Fermi parameter pp is strong, 
which allows us to perform a detailed analysis. 

Introducing in the latter model another x-variable as the ratio x = W/Mb, we may 
calculate the appropriate distribution function w{x) of the b quark in the restframe of the 
B meson as a function of the relative mass x. In Fig. 2b we plot the function w{x) for 
nisp = 0.15 GeV and various values of pp, this value for the spectator mass being frequently 
used in studies of the semileptonic decays. Note that the mass distribution shows sizable 
modihcations due to the variation of pp. Both models have the advantage of avoiding the 
mass of the heavy quark as an independent parameter. As a consequence, the phase space 
is treated correctly because they use the mesonic degrees of freedom. 

Within the approximations discussed in Section |l|, the shape of the momentum spectrum 
in the decay B —>• J/iIj + X is determined by the value of pp. As argued in Ref. |^, the 
Fermi momentum parameter is not a truly free parameter, but it is directly related to the 
average kinetic energy of the heavy quark, (p^) = through the Gaussian form of the 

Thus the value of pp can also be deduced theoretically 


distribution function of Eq. 
from a study of the b quark’s average kinetic energy. Hwang et al. calculated the 
latter in the relativistic quark model, from which they obtained pp ~ 0.5 — 0.6 GeV, this 
value being in good agreement with the one deduced from the QGD sum rule analysis of 
Ref. 1^, Pp = 0.58 ±0.06 GeV. Both results respect the limit pp > 0.49 GeV which follows 
from an inequality between the expectation value of the b quark’s kinetic energy and that 
of the chromomagnetic operator derived by Bigi et al. using the methods of the Heavy 
Quark Effective Theory 0. The lower bound from the latter inequality, however, could 
be considerably weakened by higher order perturbative corrections |0 . 

The numerical range of the Fermi motion parameter deduced from the heavy quark’s 
kinetic energy is in good agreement with a recent study of the semileptonic B decays: 
Using dilepton data, the GLEO collaboration has performed a model independent analysis 
of the lepton spectrum over essentially the full momentum range |0]. By comparing the 
results with the theoretical prediction of the AGGMM model, Hwang et al. obtained 
Pp = 0.54GeV. 
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Allowing the Fermi momentum to float within the range pp ~ 0.5 — 0.6 GeV, we 
investigate to what extent the b quark motion in the B meson can account for the observed 
spectrum in B —^ J/xfj+X decays. One might note that, in contrast to the energy spectrum 
in the semileptonic decays, the shape of the J/ifj momentum spectrum is sensitive to the 
Fermi motion parameter over a wide range of the phase space. 

Employing the spectator distribution function of Eq. (|^) , we calculated the momentum 
spectrum for the decay of B mesons produced at the T(4S') resonance. Figs. 5 and 6 show 


the comparison with the CLEO data for the set of parameters given in Section |3.2| and for 
rrisp = 0.15 GeV. In Fig. 5 one can see that (considering the fact that the AGGMM model 
does not include hadronization effects and, consequently, yields an averaged spectrum) the 
agreement of the theoretical spectrum with the data is good provided we choose the value 
Pp = 0(0.57 GeV). Again we point out the dominance of the color octet charmonium 
production mechanism by plotting also the spectrum which arises from the color singlet 
model, the latter underestimating the data by roughly a factor of three. 

To demonstrate the high sensitivity of the spectrum to the Fermi motion parameter, 
we present in Fig. 6 the spectrum for pp = 0.3 GeV (with the remaining parameters kept 
unchanged), the latter value being commonly used by experimentalists in the analysis of 
the semileptonic decays. Note the large discrepancy between theoretical prediction and 
data. Thus it is obvious that the shape of the measured J/ijj momentum distribution 
cannot be reproduced in the model, within the approximations we made, when using values 
significantly smaller than pp = 0.5 GeV. Especially the sizable contribution in the low 
momentum range requires a soft probability distribution of the heavy quark momentum. 
Further we want to emphasize that (using pp = 0.57 GeV) not only the shape but, due 
to the inclusion of color octet intermediate states in the charmonium production, also the 
normalization of the theoretical spectrum is in good agreement with the data. This way our 
analysis supports the numerical values of the NRQGD matrix elements given in Eq. (p^). 


Note that the authors of Ref. [jT^ could not account for the measured value of the 
branching ratio B{B —> J/'ip + X), because they restricted the NRQGD analysis to positive 
values of the matrix element Pq)), and because (adopting rrib = 5.3 GeV) they did 

not take into account corrections due to the bound state structure of the B meson. 

We conclude that, using the AGGMM model, the Fermi motion of the b quark can 
account for the observed momentum distribution in the B —>■ J/'ip + X decay accurately. 
The required value of the Fermi momentum pp = 0(0.57 GeV) comes out to agree with 
the range deduced from the studies of the heavy quark’s kinetic energy, as well as, with 
the range recently obtained from the semileptonic B decays. 
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We close this section with a comment on the distinct secondary bump which appears in 
the low momentum region of the measured J/ip spectrum. While this article was being 
written, Brodsky and Navarra suggested an interpretation of the latter in terms of the 
three-body decay B —>■ J/'ipKp^ in which they assumed that the underlying dynamics of the 
exclusive mode is reflected in the measured shape of the inclusive momentum distribution. 
The authors argued that the three-body hnal state accounts for a distinct enhancement in 
the inclusive spectrum, in which the position of the maximum (|k^| ~ 0.5 GeV) coincides 
with the measured bump. Thus our inclusive approach to the J/p) momentum distribution, 
when combined with the analysis of the exclusive three-body mode J/ipAp [Q, yields a 
spectrum which is in good agreement with the measurement over the whole range of phase 
space. 


5 Summary 


The analysis of the decay B J/pj + X requires the color singlet and octet contributions 
to the intermediate cc states. The octet conhgurations were found to be necessary for 


improving the predictions for the branching ratio [^|. In this article we include the octet 
conhgurations in the analysis of the J/pj momentum spectrum. We adopt the leading 
order expansion of NRQCD, where the kinematics of soft gluon radiation associated with 
the nonperturbative transitions cc —> J/pj + X is neglected. Thus we address the question 
whether and to what extent the motion of the b quark in the bound state can account for 
the observed momentum spectrum. To this end, we review in Eqs. ( 0 )-(|) and (0 the 
changes introduced to the J/pj tensor through the leading color octet conhgurations. The 
structure constants of the tensor depend on the NRQCD matrix elements {pi^'^[n\), which 
parameterize the conversion of a cc pair in angular momentum [n] and color a to the J/pJ. 
Values for the matrix elements are taken from other experiments. With these results we 
develop a formalism for the inclusive decay within the parton and the ACCMM model, 
respectively. The shapes of the spectra originating from various cc(^^^^Lj)a conhgurations 
are very similar to each other. They are sensitive to the motion of the b quark in the 
B meson. Explicit formulas for the decay spectra are given in Eqs. (IlD and (13). They 
depend on the momentum distribution function of the b quark f{x) and the Fermi motion 
p>{p) of the spectator quark, respectively. 

Numerical calculations in the parton model give a satisfactory presentation of the data, 
provided that the heavy quark momentum distribution is taken to be soft. To be explicit, 
we hnd Sp = 0(0.008 — 0.012) for the parameter of the Peterson et al. distribution function. 
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The ACCMM model can account for the observed momentum spectrum more accurately. 
The preferred Fermi momentum pp = 0(0.57 GeV) is in good agreement with the range 
Pp ~ 0.5 — 0.6 GeV deduced from theoretical studies of the heavy quark’s kinetic energy. 

Our results are obtained in the approximation of neglecting the gluon momenta at the 
hadronization stage, i.e., we include only the leading terms in the NRQGD expansion. Very 
recently, it was shown how the resummation of higher order nonperturbative corrections 
to the quarkonium production leads to the introduction of universal distribution functions 
23| . It will be a challenge for future studies to investigate the importance of these shape 
functions for the momentum distribution in the decay B —^ J/ip + X and to compare the 
theory with more accurate data. 
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Figure Captions 


Fig. 1 Separation of the various distances contributing to the B —> J/xjj + X decays. Qi ,2 
are the four-quark operators which induce the effective b —>■ ccqf transitions. 

Fig. 2a Momentum distribution function f{x) in the PM for various values of Ep. 

Fig. 2b Mass distribution function w{x) in the ACCMM model for various values of pp. 


Fig. 3 Theoretical momentum spectrum in the PM for direct inclusive J/'if) production from 
B decays at the T (dS") resonance, shown for various values of Sp and compared with 
the CLEO data; the parameters are given in Section 


Fig. 4 Same as in Fig. 3, now for Ep = 0.012 and a modihed value of the (^Pi) matrix 
element (dashed line). The dotted line shows the contribution from color singlet 
(^5i) intermediate states. 


Fig. 5 Theoretical momentum spectrum in the ACCMM model for inclusive B decays to 
J/'ijj with Pp = 0.57 GeV, compared with the CLEO data; the parameters are given 
The dashed line shows the contribution from color singlet (^S'l) 


in Section 3.2 


intermediate states. 


Fig. 6 Same as in Fig. 5, now the theoretical momentum spectrum in the ACCMM model 
for various values of pp as shown. 
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